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TEMPERATURE OF A GRAY BODY 
MOST CLOSELY FITTING 
THE SOLAR EXTRATERRESTRIAL SPECTRUM 

by 

William B. Fussell and John B. Schutt 
Goddard Space Flight Center 


SUMMARY 

A definition of "closeness" of fit of a gray body curve to the 
sun T s extraterrestrial spectrum is presented. F. S. Johnson 1 s 
data are used for the solar spectrum. For the wavelength inter- 
val 0.3 to 2.4 microns (pertinent to thermal design), the magni- 
tude of Cj calculated for use in Planck's black body formula is 
8.6504 x 10" 17 watt-cm 2 at a temperature of 5742 ± 0.5 °K. 
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TEMPERATURE OF A GRAY BODY 
MOST CLOSELY FITTING 
THE SOLAR EXTRATERRESTRIAL SPECTRUM 


by 

William B. Fussell* and John B. Schutt 
Goddard Space Flight Center 


INTRODUCTION 

The problem of simulating the extraterrestrial solar illumination for testing the thermal design 
of satellites has recently become of wide interest. If the materials for the satellite surface have 
been selected, the problem reduces to one of matching the absorbed energy from the simulating 
source— for -each material composing the surface— to the absorbed solar energy for that material. 
However, many materials have spectral absorptivity curves which vary slowly with wavelength be- 
tween wavelength limits including most of the sun’s energy. A solar simulator which matches the 
sun’s spectrum fairly well on a finely resolved basis, and which matches the total energy exactly, 
will automatically match the solar absorptivities and the simulated absorptivities closely for these 
materials. To this end, the black body function was used in constructing a gray body curve for 
which a solar fit was obtained by parametrically varying the temperature. 

METHOD 

In establishing a closeness of fit criterion for approximating the sun’s extraterrestrial spectrum 
by that of a black body, it seems pertinent at this point to recall the geometric interpretation of the 
quantities in the black body equation which can be treated as parameters. Consequently, any crite- 
rion which is established by means of parametric variations can only be defined as ultimate when it 
is the most valid one from geometric considerations. 

From the black body definition 


JxCT; C 1# C 2 ) = Cj J\(T; C 2 ) 


e C ’ AT - 1 


( 1 ) 


where k is the wavelength in cm, T is the temperature in °K, and c x and c 2 are constants (or 
parameters). Any of the quantities T, Cj, and c 2 may be chosen as a parameter. First consider 


*Mr. Fussell is now working for the National Bureau of Standards, Washington, D- C. 
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variations in C x . It is clear that any change in the magnitude of C x shifts the spectral distribution 
up or down along the ordinate depending on the sign of the chosen change for Cj. Since the parameters 
C 2 and T are geometrically equivalent (in a reciprocal manner), only one of them need be considered; 
T was chosen because of its experimental accessibility. Variations in t affect the shape of the 
black body distribution and, of course, the wavelength at which the maximum spectral intensity occurs. 
Therefore, any criterion which is to be established for fitting the sun T s extraterrestrial spectral dis- 
tribution to Equation 1 can be determined by means of parametric variations in T and C 15 or a 
gray body fit. 

The most straightforward approach is to find the temperature at which the black body maximum 
matches the sun T s maximum, 0.475 micron. This temperature is readily found to be 6060°K. But 
this method gives a black body curve which contains a greater integrated energy than the sun’s inte- 
grated spectral distribution, since by definition the black body distribution gives an upper limit to 
the emissive ability of any body at a given temperature,, Consequently, an upper bouncj on the re- 
quired temperature has been established from this matching process. Next, an estimate of the lower 
bound can be found. To accomplish this, c 1 and t may be varied in an effort to minimize the 
residual 


R = 



c 2 ) - HjdX 


( 2 ) 


where and the vertical lines represent the sun T s spectral distribution and absolute magnitude, 
respectively. From Equation 2 the lower bound is found to be approximately 5738°K. It is apparent 
that Equation 2 alone cannot give the required closeness of fit to the desired accuracy because one 
criterion cannot be used to minimize a bivariant distribution. The final criterion, the one defined as 
geometrically ultimate and, consequently, most pertinent from the thermal design viewpoint, is the 
system: 



c 2 ) - Hjd K 


and 


(3) 


^ J\(T; Cj, c 2 )d\ = JV dx . ( 4) 

The addition of Equation 4, to establish the most meaningful criterion, allows a temperature to be 
found by using matched areas as the basis. The absolute magnitude in Equation 3 is necessary in 
order to insure that the total deviation between the two curves is minimized. From Equation 1 it is 
clear that c x can be eliminated in Equation 3 by means of Equation 4, which gives a single relation 
for the determination of T . These equations give a temperature of 5742 ± 0.5°K and a value of 
8.65040 x 10 ~ 17 watt-cm 2 for C x for the 0.3 to 2.4 micron (p) region. 
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It remains to show that Equation 2 should give results identical to Equations 3 and 4, indicating 
therefore that 5742 °K is the desired temperature. If the spectral region of interest is broken into 
intervals of J^(T; c 1 ) > and J\(T; C 1 ) < H^, where the intervals of the first kind are designated by 
Aj t and those of the second kind by 1 , then Equations 2 and 3 may be rewritten: 


f |J\(T; C,) - H,Jdk 

+ f IJ\(T; Cj) - HjdX 

k.i 

-kj.i 

= f c i J\(T)d^ - 

f d\ - f Cj J\(T)dX + f 

Ja i,i 

K.l . 


(5) 


where the first subscript of A refers to the interval type and the second is an index which enumer- 
ates subintervals. When the derivatives of the right-hand side of Equation 5 (ignore Equation 4) are 
taken with respect to c x and T in an effort to obtain simultaneous equations for C x and T, the 
following results are obtained: •- 


i 


*1 , 1 » A 2 . 1 


J\(T)d\ 


0 


and 


c i 


i 


* 2.1 


BT J\( T ) 


0 . 


Consequently, Equation 4 must be included in order to define c x and obtain the desired results. This 
implies the equivalence of Equation 2 and the system of Equations 3 and 4. By taking the derivative 
of the right-hand side of Equation 5 with respect to T and equating the result to zero in order to ob- 
tain a minimum (make use of Equation 4), Equation 5 may be rewritten: 

|J\(T; c,) - Hpjd K + f I J\(T; c,) - HjJd X < £ , (6 ) 

J *i,i Ja 2,i 

where e is a number characterizing the fit of the black body curve to the sun' s distribution and 
Aj , and A 2 j become the sets {Aj_ .} and {a 2 .]■ , respectively, with l < i < N and l < j < M. 

If the fit is exact, e = o and the equality sign holds; however, if the fit is not exact, the left-hand 
side of Equation 5 gives the best fit, which means that Equations 3 and 4 give the desired temperature. 

RESULTS 

The results given after each method of calculation described in the previous section are for the 
wavelength range of 0.3-2. 4 jn. Calculations were also made for the intervals 0. 4-1.2 and 0.220-7.0 u ; 
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the former region is of interest for purposes of solar cell illumination, and the latter is a complete 
cover of Johnson's data.* All results are given in Table 1. The auxiliary range of 0.3-2. Op is in- 
cluded for comparison purposes. 

It is apparent that the best fit (2.97 percent) was obtained for the 0.4-1.2p region, which indi- 
cates that the adjusted black body curve conforms reasonably well to Johnson's data throughout the 
visible and into the near infrared. Less satisfactory agreement (6.08 percent) was obtained for the 
region of 0.3-2. Op because a single black body formula cannot be made to fit the steep rise in emis- 
siveness characteristic of the sun's radiation in the region of 0.34-0.4p. Agreement is enhanced, 
however, when the region is extended to 2.4p because the additional 2.0-2.4p region, where the fit 
is good, makes the subvisible region (0.34-0.4p) of somewhat less weight. The fit for the complete 
range covered by Johnson's data, 0.220-7. Op , is least satisfactory. The black body curve (with only 
two adjustable parameters) cannot be made to fit the steep rise in the 0. 220-0. 4p region as well as 
the tail of the sun’s spectral distribution in the infrared. 

In order to clarify the fit to the solar extraterrestrial spectrum by the above gray body matching 
technique, a gray body curve was calculated for the 0.3-2.4p region. The gray body curve of this 
region for a temperature of 5742°K is compared with the solar data in Figure 1. According to col- 
umn 6 of Table 1 the areal mismatch of these curves is 5.94 percent when computed relative to the 
solar constant. It can be concluded that a gray body curve agrees well with the solar data curve, 
except for the irregular region of 0.41-0.55p. What has happened is that virtually exact fits occurred 
at 0.43 and 0.52 p, and the intermediate irregularities could not be fit by a single gray body curve. 

Table 1 


Results of Gray Body Fit to Solar Extraterrestrial Spectrum. 


Range 

(microns) 

T (°K) 

C x (watt-cm 2 ) 

/ Hx.d\ (watt/cm 2 ) 

J l J \ " Hjd\(watt/cm 2 ) 

' |J ^ |d \,00 
/ H^d K 

0.4 - 1.2 

6142 ± 0.5 

6.76111 x 1CT 17 

9.86137 x 10" 2 

2.92699 x 10'3 

2.97 

0.3 “2.0 

5730 ± 0.5 

8.72019 x 10‘ 17 

1.29419 x 10" 1 

7.86669 x 10~ 3 

6.08 

0.3 -2.4 

5742 ± 0.5 

8.65040 x 10- 17 

1.32629 x 10" 1 

7.88560 x 10~ 3 

5.94 

0.220 “7.0 

5693 ± 0.5 

8.80334 x 10~ 17 

1 .39424 x 1 0 _1 

1.03376 x IQ" 3 

7.41 


DISCUSSION 

If the black body approximation to the solar spectrum is to be useful from the thermal design 
viewpoint, the average absorptance of a spacecraft coating material, as calculated from the black 
body spectrum, should be fairly close in value to the average solar absorptance of the material. 
Mathematically expressed, it is desirable that a BB - a s , where 


•Johnson, F. S., "The Solar Constant," /. Meteorology 11(6): 431-439, December 1954. 
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Figure 1 -Comparison of a 5742°K gray body curve with the solar extraterrestrial spectru 




J o 


J\(T)d\ 


/•CD 

a\H\ dA. 



with 


a BB = the average absorptance for a black body, 
a g = the average solar absorptance, 

H\ = the solar spectral intensity per unit wavelength bandwidth, as given by Johnson, 
J\(T) = the Planck black body spectral intensity for absolute temperture T, 

= the absorptance of the material at wavelength K . 


For simplicity, it is desirable to normalize and to H\ and 


= 


/»C0 

Jo 




Jk 




•l 


h & 


Equations 7 and 8 now become 


j 'OD 

k 

0 


\Jx.(T)dX 


( 9 ) 


and 


/» 00 

Jo 


a ^d A. 


( 10 ) 


Arbitrarily, a 10 percent proportional difference in a BB and a s will be chosen as the criterion 
of usefulness; that is, if 


I a BB 3 s I < 

V s E 1 0.1 , 

p ’ 


(ii) 


the black body approximation will be considered useful 0 Equation 11 is actually 


/•CD 

j 0 [j^(T)dX - H^]a x d\ 


r 

Jo 


d\ 


* 0.1 


( 12 ) 
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It can be seen that the results given in column 5, Table 1 for the quantity 


I 


IjjtfT) - H£|dX - A p 


( 13 ) 


do not indicate the magnitude of E p (Equation 11). Consider now the following two cases: 
1. If a x is approximately constant, then the proportional error E p becomes 


a X I £ [Km - dxl 


rco 

H x 

Jo 


(14) 




and, since 


J 'OO /'CO 

Jx(T)dX = H^dX 
0 Jo 


(15) 


it is clear that 


i IA< T > - "i] 


dk = 0 


(16) 


and E = o . 

p 


2. If can be represented by a power series in k , that is, 


da x\ (X-X 0 ) 2 /d 2 a x 


dk 2 


*k ' ( a \) + i k ~ X o)("d)T n + 2! 


(17) 


then 


1 1 k( T > - Hx] 

poo 

H£axdX 

Jo 


(18) 


in which the error is now weighted by the slope. 
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"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 
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